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Variations in the apparent flux and position of the radio source Cygnus-A were recorded 
at 53 and 108 megacycles per second using a two-element, phase-sweeping interferometer 
located at Boulder, Colo. An ionospheric sounder operating at Ellsworth, Nebraska, 
provided, for a few hours each day, simultaneous vertical-incidence measurements on the 
ionosphere at its intersection with the line of sight from Boulder to 1 he radio s1 ar. Amplitude 
scintillations observed at Boulder over a twelve-month period are compared with ionograms 
taken at Ellsworth. Positive correlation is found between amplitude scintillations and 
spread F, while no significant correlation is found with sporadic E. 

Detailed analysis of the scintillations indicates thai the probability distribution of the 
amplitude can be represented by the Rice probability distribution function. The zenith- 
angle dependence of the amplitude scintillations does not agree with a theory based upon 
isotropic ionospheric inhomogeneities. 

The root-mean-square value of angular scintillations is proportional to the square of the 
wavelength, in accord with a theory of diffraction by ionospheric irregularities. Comparison 
of angular scintillations with amplitude scintillations indicates that, for elevation angles of 
15° to 50°, the region of the ionosphere responsible for scintillation lies near- the border be- 
tween the Fresnel and Fraunhofer diffraction regions for both frequencies. 

Slow, irregular angular- variations are commonly observed in the daytime at both 
frequencies. These variations are attributed to lens-like ionospheric irregularities having 
dimensions as large as 200 kilometers. 



1. Introduction 

It has been recognized since 1950 that the apparent 
fluctuations in intensity of discrete radio sources 
arise in the ionosphere. In a manner analogous to 
the twinkling mid dancing of visible stars, these 
"amplitude scintillations" are accompanied by "an- 
gular scintillations", which are fluctuations in the 
apparent position of the radio source. 

Radio-star scintillations offer a convenient means 
for studying the irregular structure of the ionosphere. 
The F region is generally believed to play an impor- 
tant role in the production of scintillations, but our 
knowledge of the effect of the E'regio i is confused by 
incomplete and apparently conflicting evidence. One 
of the purposes of the experiment reported here was 
to improve our understanding of the ionospheric 
conditions and levels responsible for scintillations. 
The experiment was designed to be particularly 
sensitive to any relationship between scintillations 
and the E region. The results have indicated little 
or no such relationship, but do show a connection 
between scintillations and spread echoes in the F 
region. In addition to the comparison of scintilla- 
tions and ionospheric phenomena, detailed analysis 
has been made of the scintillations themselves. In 
connection with this, some statistical aspects of 
amplitude scintillations are found to disagree with a 
theory based upon plane-earth geometry and spher- 
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teal irregularities whereas the statistical variations 
in angular scintillations agree with this theory. 

A third aspect of the experiment has been the 
observation and study of slow angular variations in 
the apparent position of the source. These varia- 
tions are due to refraction by individual large irregu- 
larities, in contrast to amplitude and angular scintil- 
lations which are due to the diffraction effects of a 
number of smaller irregularities. 

A description of the planning and instrumentation 
for the experiment, along with a brief summary of the 
pertinent literature, has already appeared [Lawrence, 
1958]. In addition, two review papers [Booker, 1958; 
Little et ah, 1956] are available which outline present 
knowledge concerning scintillations and which de- 
scribe the usefulness of scintillation observations as a 
tool for studying the irregular structure of the 
ionosphere. 

2. Description of the Observations 

The radio source Cygnus-A has been observed at 
Boulder (latitude 40°7'.7 N, longitude 105°14'.3 W) 
at frequencies of 53 and 108 Mc/s, using a two-ele- 
ment, phase-sweeping interferometer. The phase- 
sweeping system permits simultaneous observation of 
the apparent flux from the radio source and of the 
relative phase of the radio waves at the two antennas. 
Amplitude scintillations are simply variations in the 
flux, while angular scintillations may be deduced from 
the known spacing of the antennas, the general direc- 
tion of the radio source, and the variations in relative 
phase of the radio waves. 
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From February 1958 through February 1959 ob- 
servations were obtained regularly during the three or 
four hours each day following the rising of Cygnus-A. 
Figure 1 shows how the time of such observations 
varies throughout the year and indicates the seasons 
during which the effective ionosphere was sunlit. 

The two antennas are separated by 475 m on a 
horizontal, east-west baseline and the resulting inter- 
ferometer lobes are of nearly constant spacing over 
the region of the sky in which the observations were 
made. At a radiofrequency of 108 Mc/s the lobe 
spacing is about 0.4 deg; at 53 Mc/s it is twice as 
great. The measured component of angular scintil- 
lations is, of course, perpendicular to the lobes. This 
means that, near the horizon, angular scintillations 
can be observed only in the horizontal direction, 
while at an elevation angle of about 35° the horizontal 
and vertical components of angular scintillations are 
observed with equal sensitivity. 

At the time when Cygnus-A is about 15° above 
the north-east horizon at Boulder, the line of sight 
intersects the E region over Ellsworth, Nebraska. 
During the appropriate portion of each day, an iono- 
spheric sounder was operated at Ellsworth at 5-min 
intervals. As the radio source rises, the intersection 
point of the line of sight with the E region moves to- 
ward Boulder, and the intersection with the F region 
passes within about 100 km of Ellsworth. Figure 2 
is a map of the Boulder-Ellsworth area showing the 
plan view of the intersection point for various iono- 
spheric heights. The numbers along the curves indi- 
cate the local sidereal time (LST) at Boulder. 



3. Instrumentation 

The instrumentation for the scintillation measure- 
ments has already been described in some detail 
[Lawrence, 1958], but is briefly reviewed here. The 
two frequencies, 53 and 108 Mc/s, are observed with 
the same antennas and with similar receiving and 
recording systems. 

Antennas. The two 40-ft paraboloidial reflecting 
antennas are modifications of standard 28-f t parabo- 
loids manufactured by the D. S. Kennedy Company. 
The 6-ft radial extensions were designed and con- 
structed at NBS. The antennas, separated b} r 475 
m, are equatorially mounted and equipped with 
sidereal drives for tracking radio stars. The dipole 
feed, with four reflecting elements, is designed for use 
simultaneously at both frequencies. A multicoupler 
permits the use of a single coaxial transmission line 
from each antenna to the central receiving shelter. 
No preamplifiers are used at the antennas. 

Receivers. At each frequency, the dual receiving 
channels are connected together as a Ryle inter- 
ferometer [Ryle, 1952]. A rotary phase shifter has 
been added to one channel so that the interferometer 
pattern may be swept across the sky. This phase 
shifter, usually driven at a uniform rate of 10 rpm, 
converts the system to a "phase-sweeping" inter- 
ferometer. A "phase-tracking" feature was in- 
cluded in the original design, whereby the inter- 
ferometer pattern was servo-controlled to follow the 
apparent motion of the radio source. This feature 
permits the direct recording of angular scintillations. 
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Figuke 1. Time of observation of Cygnus-A, and state of illumination of the ionosphere. 
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Figure 2. The variation, with local sidereal time at Boulder, 
of the position at which the line-of -sight from Boulder to 
Cygnus-A penetrates the ionosphere. 

The three curves apply to ionospheric heights of 90, 100, and 300 km. 

However, it keeps a zero-crossing of the inter- 
ferometer pattern always on the source and so 
prevents recording of amplitude information. As 
a result, it is difficult to distinguish authentic angular 
scintillations from the random behavior of the servo 
loop which accompanies deep amplitude fades. 
Accordingly, all of the observations reported here 
were made in the phase-sweeping mode. 

The receivers have bandwidths of 500 kc/s, and 
were operated with output time constants of one 
second. The double-conversion superheterodyne 
system uses intermediate frequencies of 31.5 and 
7.5 Mc/s. 

Digital recorder. In order to permit the use of a 
computer for processing the data,, chart recordings 
were supplemented by simultaneous digital record- 
ings on punched paper tape. A simple, meter-type 
analog-to-digital converter sampled the output from 
each interferometer once per second, and punched 
the two outputs alternately on a single paper tape. 
The six-bit binary code provides a full-scale accuracy 
of 1.5 percent, approximately the same as careful 
manual scaling of a pen recording. 

4. Data Analysis 

The following is a summary of data-analysis 
procedures. The appendices and the following 
paper by the same authors, contain a detailed ex- 
planation of the procedures outlined in this section. 

The observations made in the progress of this 
experiment consist of scintillation records from the 
phase-switch, phase-sweep interferometer located at 
Boulder, and of ionograms from an ionospheric 
sounder located at Ellsworth, Nebr. 

Ionograms. Vertical-incidence, sweep-frequency 
ionospheric records were taken at 5 min intervals 



for 1 hr, 1230 through 1330 LST each day during 
the time that the ray path from Cygnus-A passed 
through the E region (100 km) near Ellsworth. 
After this period, ionograms were made each half- 
hour for the next 3 hr. These ionograms have been 
scaled for critical frequencies, minimum virtual 
heights, and for the occurrence of special features, 
particularly for the absence or presence of spread F 
and of the various types of sporadic E. Each type 
of information from the ionogram was converted to 
an appropriate single-digit scale and punched in a 
separate column of a uniquely identified IBM card. 
Simultaneous scintillation information, described 
below, was later added to each card. 

Scintillation records. The scintillation records, 
both analog and digital, consist of a sine wave ("car- 
rier") having a frequency determined b} r the rate 
of lobe sweeping and an amplitude and phase which 
are modulated by the ionosphere. For the present 
experiment several methods were used to extract, 
from the digital records, the amplitude and phase 
modulation. However, all of these methods depend 
upon the fact that the envelope of the carrier is 
affected only by ionospheric amplitude effects and 
that the times of zero -crossings are affected only 
by ionospheric phase effects. Intermediate points 
on the carrier are affected by both. The power 
index 
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have been computed to measure amplitude effects. 

P is the power from the discrete source and may be 
read directly from the envelope of the carrier since a 
square-law detector was used in the scintillation 
receivers. A = ^P, and a bar denotes an average 
over the record. The reason for computing both 
of these quantities will be shown in a later section. 

The variance of the angular scintillations, 2 , was 
computed as a measure_of phase effects. For any 
given elevation angle, B 2 can ho deduced from the 
observed value of a a , the standard deviation of the 
phase difference, a, between the two ray paths to 
the interferometer. Variations in a may be deter- 
mined directly from the variations in the interval 
of time between successive zero-crossings. (See 
following paper by the same authors.) 

The first attempt to measure amplitude scintilla- 
tions was to assign visually one of five possible cate- 
gories to each 5-min section of the scintillation 
records (see appendix 1). Those sections for which 
there were no visible variations in the envelope of 
the carrier were assigned a "depth index' 7 of 1, those 
for which there were shallow variations in the enve- 
lope were assigned a depth index of 2, etc., to the 
depth index 5 for records with amplitude effects so 
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deep that the lobes disappeared. Several examples 
of each of these categories of records were scaled 
by hand so that an estimate could be made of the 
range in values of (AA/A) 2 and of (AP/P) associated 
with each category (see the table in appendix 1). 

A second approach was to recover P and a from 
the carrier and then to plot these two quantities as 
functions of time. The resulting graphs are referred 
to as power plots and phase plots respectively. 
Two different methods were used ("zero-crossing" 
and "filter") to make this recovery (see following 
paper). The power and phase plots may be used to 
compute the quantities (AP/P), (AA/A) 2 , and 8 2 , 
and it is also possible to use them to determine the 
power spectra of amplitude and angular scintillations. 

Some of the records, because of large amplitude 
and/or phase effects, are not amenable to analysis 
by either the "zero-crossing" or "filter" methods. 
These records were analyzed using a statistical ap- 
proach which does not give detailed power and phase 
plots but which gives the quantities P n and <r d , 
quantities which may be related to (AA/A) 2 and 6 2 
(see following paper by the same authors). 

5. Results and Discussion 

The relation between amplitude scintillations and 
ionospheric parameters. Figure 3 displays the distri- 
bution with elevation angle of the Ellsworth iono- 
grams and the coincident scintillation records which 
were used to investigate the correlation between the 
two types of observations. Figure 4 shows a definite 
correlation at both radiofrequencies between the 
depth of amplitude scintillation and the occurrence 
of spread F. Here the Fisher Z r function [Walker 
and Lev, 1953] is plotted at the weighted mean 
elevation angles of the 10-deg class intervals shown 
in figure 3. The upper scale gives, for each eleva- 
tion angle, the ionospheric height at which the line 
of sight to Cygnus-A comes closest to the zenith at 
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Figuee 3. The grouping for statistical analysis, with elevation 
angle and sidereal time, of the scintillation data. 

The "effective ionospheric height" is the height at which the line-of-sight 
from Boulder to Cygnus-A comes closest to the zenith at Ellsworth, Nebr. 



Ellsworth. N is the number of samples used in the 
computation of each value of Z r , and the vertical 
bar through each point has length 2o- Zr , twice the 
expected standard deviation of the correlation co- 
efficient. Each graph is derived using the total data 
sample. Division of the data into sunlit and dark 
ionospheric groups results in correlations (not shown) 
which have, in general, a value of Z r which is 0.1 
higher (dark ionosphere) and 0.1 lower (sunlit iono- 
sphere) than the correlations shown, with corres- 
pondingly larger values of o- 2r . 
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Figure 4. The correlation between amplitude scintillations 
and spread F, as a function of source elevation for each 
radio frequency. 

N is the number of cases. 



336 



In the same manner, figure 5 displays the correla- 
tion, using the total data sample, between the depth 
of amplitude scintillations and the occurrence of 
sporadic E. Division of the data according to sunlit 
and dark ionosphere, as well as by season, resulted 
in no significant change from the results shown. 

Other ionospheric parameters showed little or no 
correlation with the depth of amplitude scintilla- 
tions. In particular, there was no significant corre- 
lation, at either racliofrequency, with joF2, joEs, 
h'Es, or / min . Correlation with h/F, the minimum 
virtual height of the F region, resulted in a Z r of 
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Figure 5. The correlation between amplitude scintillations 
and sporadic E, as a function of source elevation for each 
radio frequency. 



about 0.2 to 0.3 at 53 Mc/s and a Z r of about 0.1 to 
0.2 at 108 Mc/s, at i^-region effective heights. 

Many, but not all, previous observers have found 
scintillations to be correlated with the phenomenon 
of spread-/' 7 ionospheric echoes. In the early days 
of scintillation study, Ryle and Hewish [1950] and 
Little and Maxwell [1951] both found such relation- 
ships, while Mills and Thomas [1951] found positive 
correlation with an "jF-region index". Later Bolton, 
Slee, and Stanley [1953] reported positive correla- 
tion during the winter months; as did Wild and 
Roberts [1956a] for nighttime scintillations, although 
in both cases the ionosonde was located many 
hundreds of kilometers from the point of i^-region 
penetration. Both Duefio [1955] and Hartz [1958] 
failed to find correlation with ionosonde data taken 
at points appreciably far from the optimum for F- 
region penetration, although Hartz and, more re- 
cently, Dagg [1957] noted that the diurnal variations 
are similar. Briggs [1958], on the other hand, has 
recently studied the correlation between scintilla- 
tions and spread F to deduce the spatial extent and 
the height of the region responsible for scintillations. 
Wright, Koster, and Skinner [1956]; and Koster, 
[1958] have noticed a quite high positive correlation 
with equatorial spread F, but this may be a differ- 
ent phenomenon from the spread F at temperate 
lat itudes. 

The diurnal variations of spread F and scintilla- 
tions are generally reported to be similar, and a 
comparison between figures 6 and 7 indicates that 
this is clearly the case for the present experiment. 
This fact alone might explain some of the correla- 
tions which have been reported. In other cases, a 
disadvantageous location of the ionosonde may have 
reduced or destroyed the correlation. We have at- 
tempted to avoid both pitfalls, first, by removing 
the average diurnal and seasonal variation from both 
the scintillation data, and the spread F observations 
(see appendix 2), and second, by locating the iono- 
sonde at Ellsworth, Nebr. Although Ellsworth was 
chosen specifically for E'-region observations, the 
location should provide quite reliable /''-region data 
for elevation angles of Cygnus-A in the range from 
25 to 65 deg above the northeastern horizon at 
Boulder. 

The positive correlation between spread F and 
the depth of amplitude scintillations, as shown by 
figure 4, agrees with many of the results described 
previously. The location of the ionosonde, although 
not ideal for F-region effects, is quite satisfactory if, 
as has been reported [Briggs, 1958; Koster, 1958; 
Yeh and Swenson, 1959; and Reid, 1957], the irregu- 
larities responsible for spread F occur in patches 300 
to 500 km in lateral extent. The data revision proc- 
ess, which removed the diurnal and seasonal varia- 
tion from the ionospheric and scintillation param- 
eters, insures that any correlation which remains 
must be the result of detailed hour-to-hour or day- 
to-da}^ correlation between the two variables (see 
appendix 2). 

The correlation, for a given elevation angle of 
Cygnus-A was generally higher at 53 Mc/s than at 
108 Mc/s. This result might have been anticipated 
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Figure 6. The frequency of occurrence, averaged over a year, 
of spread F and of sporadic E at the times of scintillation 
observations. 

since the refractivity of the ionosphere varies in- 
versely with the square of the radio-wave frequency. 
The nearly constant value of correlation, at either 
radiofrequency, over the whole range of elevation 
angles, indicates that the patches of spread F cover 
large areas. This result has the unfortunate effect 
of making height determination, by the method used 
in this experiment, impractical. The correlation 
was considerably higher at nighttime than in the 
daytime, and this result tends to reinforce the idea 
that scintillations are related to spread F, since 
spread F reaches a maximum shortly after local 
midnight. It is possible that this variation of cor- 
relation with local time is partly a result of our 
inability to remove completely the effect of diurnal 
and seasonal trends. It is clear that correlations of 
this type could be found simply by comparing trends, 
even though there is little, or no, detailed correlation 
present. 

The correlation between scintillations and spread 
F is significant, but still far from perfect; however, 
two features of our experimental technique have 
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Figure 7. Hourly average values of scintillation depth index 
and of the statistical quantity P n , both converted to (A A I A) 2 
and to the "noise-to-signal ratio", k. 

tended to reduce the apparent correlation. First, 
Ellsworth is not ideally located for i^-region com- 
parisons; second, and perhaps more important, we 
have used a spread-7* 1 index which depends only 
upon the absence or presence of spread F, not upon 
its intensity. By themselves, these two effects 
would be insufficient to reduce an otherwise perfect 
correlation to the observed value of about 0.3. 
Therefore, it seems probable that scintillations 
sometimes arise in ionospheric regions, such as E or 
upper F, which do not produce observable spread-.? 7 
echoes. 

A number of scintillation observers have searched 
for correlation with E'-region phenomena. Mills 
and Thomas [1951] found no correlation with 
sporadic E, whereas Dueno [1955] and Bolton, Slee, 
and Stanley [1953] reported a general, not detailed, 
positive correlation between the occurrence of 
sporadic E and the occurrence of scintillations. In 
each case, the ionosonde was located more than 200 
km from the ^-region penetration point. Wild and 
Roberts [1956a] found positive correlation with day- 
time fEs, but no correlation at night. Hartz [1958] 
working in Canada, found a general negative corre- 
lation with sporadic E. 

Because sporadic E seems to occur in smaller 
patches than does spread F, it is significant that the 
Ellsworth station was located at precisely the correct 
place to study correlations at ^-region heights. 
The ray path from Cygnus-A to Boulder passes very 
near the zenith at Ellsworth at the 100 km level and 
is within 100 km of the zenith for elevation angles in 
the range from 10 to 22 deg above the northeastern 
horizon at Boulder. As mentioned previously, 
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the data revision process retains any correlation 
between detailed events but markedly reduces the 
effect of diurnal and seasonal trends. The results 
of other workers [Hartz, 1958; Dyce, 1955] as well 
as those shown in figures 6 and 7 indicate that the 
diurnal variation of the occurrence of sporadic E 
and the diurnal variation of the occurrence of scintil- 
lations are inversely correlated. 

Interim results [Lawrence and Jespersen, 1959] of 
the present experiment indicated a negative correla- 
tion (at i^-region heights) between scintillations and 
sporadic E. Those early results were obtained from 
only the first six months (Feb. 58-July 58) data, 
and the diurnal and seasonal trends had not been 
removed. It appears in retrospect that these corre- 
lations were produced by the inverse diurnal trends 
in the two variables. They, as well as the present 
results, are derived from data taken at 5-min 
intervals in the relevant range of elevation angles. 
Accordingly, the sensitivity of the result to any real 
relationship should be high. 

The present results, shown in figure 5, indicate no 
significant correlation, at either radiofrequency, with 
sporadic E. Although this contradicts the weak 
negative correlation mentioned in the interim report, 
we infer that there is no observable, detailed con- 
nection between sporadic E and amplitude scintilla- 
tions at this latitude. 

Excepting spread F\ the only ionospheric param- 
eter to exhibit significant correlation with scintilla- 
tions was h'F, the minimum virtual height of the 
/'region. This positive correlation was higher at 
53 Mc/s than at 10S Mc/s and tended, even more 
than did the correlation with spread F } to be re- 
stricted to a range of elevation angles which indicated 
/^-region heights. This suggests that the properties 
of the electron-density profile which determine h'F 
are related to scintillations, but do not have good 
serial correlation over such large distances as does 
spread F. In addition to this, the correlation may 
be related to the diurnal and seasonal variation of 
radio star scintillations, a variation which has been 
the object of much study [Bolton, Slee, and Stanley, 
1953; Hartz, 1958; Dagg, 1957; Briggs, 1958; 
Koster, 1958; Yeh and Swenson, 1959]. Figure 7 
shows the observed variation in (AA/A) 2 . This 
variation contains a mixture of both seasonal and 
diurnal effects since only one source, Cygnus-A, was 
used. Theory indicates [Booker, 1958] that, in the 
Fresncl (near zone) region, (AA/A) 2 is proportional 
to the square of the distance, z, between the observer 
and the diffracting screen (see eq (2)). If h'F is 
related to z, then one would expect a correlation 
between amplitude scintillations and h'F; also, any 
seasonal or diurnal variation in h'F would produce 
a cor re ponding variation in (AA/A) 2 . However, 
hmaxF is probably a better measure of z than is 
h'F. In view of this possibility, the observed diurnal 
and seasonal variation in h mSiX F, at this latitude for 
the period corresponding to the times that scintilla- 
tion data were taken, was used to obtain zforeq (2), 
and the resulting variation in (AA/A) 2 was com- 
pared with the observed variation (fig. 7). At 53 
Mc/s (using the P n curve) the agreement was quite 



good. No comparison was made at 108 Mc/s since 
the available portion of the P n curve shows little 
variation. The agreement at 53 Mc/s shows that 
a hypothetical diffracting screen which always re- 
mains at the height of maximum electron density, 
and whose statistical properties never change, will 
produce the same diurnal variation in amplitude 
scintillations as is observed at Boulder. However, 
such a model seems quite unreasonable since it 
implies that all the irregularities are forced to move 
with h max F, and that their intensity, AN, is inde- 
pendent of electron density, N. 

Amplitude scintillations. The power plots ob- 
tained from the "zero-crossing" and "filter" methods 
provide a convenient basis for detailed analysis 
of amplitude scintillations. Using the method 
described by Tukey [1949], power spectra extending 
over the frequency range 0.005 to 0.17 c/s have 
been calculated from the power plots. Four such 
spectra, together with the corresponding values of 
(AA/A) 2 , are displayed in figure 8. The ordinate 
value is (he common logarithm of the variance per 
unit frequency interval, expressed in units of the 
square of the average power. Since the spectra 
come from power plots, and these cannot be ob- 
tained for the highly disturbed records, the two 
upper curves in figure 8, although they are the 
highest spectra available, should be taken as typical 
of only moderate scintillation activity. The two 
lower spectra were derived From data with little 
scintillation activity, and are close to the noise 
level. 




O.08 
FREQUENCY, 



0.12 
c/s 



Figure 8. The figure shows the highest and lowest observed 
amplitude power spectra for each radio frequency. 

The lower level at each frequency indicates the background "noise" level. 
Tlie units of IF are the variance per unit frequency interval expressed in units of 
the square of the average power. 
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Each pair of spectra in figure 8 was derived from 
simultaneous power plots at the two radiofrequen- 
cies. Comparison of these spectra, and of many 
others like them, indicates that there is no detailed 
correlation in the frequency domain. (The quasi- 
periodic variation superposed on the spectra is 
caused by the finite length of the data samples 
and should not be misconstrued as an indication 
of correlation between the spectra [Blackman and 
Tukey, 1959].) No correlation is to be expected 
for scintillations which result from a diffraction 
process because the distance, measured in wave- 
lengths, between the observer and the diffracting 
screen is very different for the two radiofrequencies. 
Refraction, on the other hand, should cause detailed 
correlation. Thus, we conclude that, at 108 Mc/s, 
refraction plays no observable role in the produc- 
tion of the amplitude scintillations we observe, 
in contrast to the observations made at lower fre- 
quencies by Wild and Roberts [1956b.]. The lack of 
correlation also indicates that there are no correlated 
instrumental errors, such as might be introduced 
by variations in the power-line voltage or by move- 
ments of the antennas. 

Excepting the short-period fluctuations which 
are the effect of finite sample length, the power 
spectra follow a definite pattern. Typically, log 
W decreases linearly with frequency, up to a fre- 
quency of about 0.1 c/s, and becomes nearly con- 
stant thereafter. The negative slope of the decreas- 
ing portion of the spectrum is directly related to the 
"fading period," t, which was used in a previous 
study of the amplitude scintillations [Lawrence and 
Jespersen, 1959]. (In that study, t was defined as 
twice the average interval between crossings of the 
average power level.) Although the upper example 
in figure 8 does not seem to concur, the negative 
slopes of pairs of spectra derived from simultaneous 
scintillation records made at the two radiofrequen- 
cies, are generally very similar. Such independence 
between fading rate and frequency is to be expected 
for scintillations which arise from a diffraction 
process. 

The power spectra are nearly flat for frequencies 
greater than 0.1 c/s, but the shape of this portion of 
the spectrum is severely influenced by the discrete 
sampling of the original data. When a function, 
x(t), is sampled at equal intervals of time, At, then 
the frequencies 2f N ±J, 4.f N ±f, 6/^±/, etc. are not 
distinguishable from the frequency / [Blackman and 
Tukey, 1959]. Here, f N , the Nyquist frequency, is 
defined as f N = ){At. Thus, if frequencies higher than 
jx are present in the data they will contribute their 
energies to lower frequencies in the spectrum. 
This phenomenon is frequently described as 
"aliasing" or "folding" of the higher frequencies into 
the range of interest of the power spectrum. 

As mentioned previously, the power spectra could 
be computed onlv for records which displayed mod- 
erate, slow scintillations. On such records, most of 
the variance due to scintillations is composed of 
frequencies less than j N) and those frequency com- 
ponents greater than j N will be largely due to white 
noise. The effect of white noise, when folded back 



into the spectrum, has been examined and is found 
to increase the general level of the spectrum by a 
nearly constant amount. In order to avoid aliasing 
the power spectra, one would ordinarily use appro- 
priate low-pass filters to remove the higher frequen- 
cies before sampling. However, in the case of the 
"zero-crossing" method (see following paper by the 
same authors), the necessary filter would have an 
upper frequency cutoff equal to the phase-sweeping 
frequency and would thus remove the signal. 

The power plots have also been analyzed in terms 
of {AAjA) 2 and (AP/P), the amplitude and power 
indices [Lawrence and Jespersen, 1959]. The rela- 
tion between (AA/A) 2 and (AP/P) depends upon the 
distribution law obeyed by the amplitude, A, of the 
radio wave. Thirty-one records selected at random 
were used to produce the dots shown in figure 9, 
where it may be noted that the observed values of 
amplitude index (AA/A) 2 and power index (AP/P) lie 
close to a smooth curve. 

Since the relation between these two quantities 
depends upon the distribution function of the 
amplitude of the radio wave, this result gives us 
some information about the distribution functions 
themselves. Any particular probability density 
function of amplitude should produce a single value of 
(AA/A) 2 and a corresponding single value of (AP/P). 
The fact that the points on figure 9 lie near a single 
smooth curve implies that the family of probability 
density functions from which our scintillation samples 
have been taken can be described by a single param- 
eter, the position along the line depending upon the 
value of that parameter. 

The solid line in the figure shows the theoretical 
relation which would occur if the observed signal 
were composed of a steady (undiffracted) component 
and a random component consisting of a number 
(three or more) of independent, randomly phased, 
diffracted waves. The probability density function 
from which this curve has been derived was first 
described by Rice [1944 and 1945] and has been 
discussed in terms of its application to ionospheric 
problems by Ratcliffe [1956]. This probability 
density function, which we shall call the "Rice 
distribution", is really a family of distribution 
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Figure 9. The relationship between amplitude and power 
indices, compared to the theoretical curve from the Rice 
distribution. 
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functions para metric in the quantity 



a= 



amplitude of the steady component 



rms instantaneous amplitude of the random component 

V2(rms instantaneous amplitude of the steady component) 
rms instantaneous amplitude of the random component 



Norton et al., [1955] have defined the quantity 



, 2 _ , 2_r rnis ^stantaneous amplitude of the random comj 
L rms instantaneous amplitude of the steady compo 



)onent~~] 2 
onent J 



which has the advantage of being a ratio of mean 
squares, k may be thought of as the rms noise-to- 
signal ratio of the radio wave. In the present 
application "signal" means the steady component 
unaffected by the irregularities of the ionosphere, 
and "noise" means the combination of the various 
diffracted waves. 

The Rice distribution degenerates into the familiar 
Rayleigb distribution in the limit as k approaches 
infinity, and approaches a displaced Gaussian dis- 
trihntion as k approaches zero. The details of the 
derivation of this curve have been given elsewhere 4 
[Lawrence and Jespersen, 1959]. 

The apparent agreement between the amplitude 
scintillation records and the Rice distribution means 
that the arbitrary amplitude and power indices 
may be replaced by the physically more meaningful 
parameter k. (This agreement is also justification 
for the method chosen to produce artificial test data, 
as described in the following paper.) Each depth 
index naturally covers a considerable range of values 
of k. Extreme and average values are approximately 
as shown below. 
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The family of Rice distributions as defined by 
Rice [1944 and 1945] applies to the amplitude of the 
radio wave and is thus quite convenient for the 
study of communications signals. For radio-star 
observations it is usual to employ a square-law de- 
tector and thus obtain records which are linear in 
power. The Rice distributions, G (P), plotted in 
terms of power, are shown in figure 10. These curves 
accurately match the observed distributions of 
chart deflections. On these curves the long tail 
which exists for large values of k is particularly 
noticeable. The corresponding large deflections of 
the record are indeed prominent, and may easily 
give the erroneous impression that they are caused 



unrelated to random 



by some special mechanism 
diffraction. 

Many simultaneous pairs of records at 53 Mc/s 

and 108 Mc/s have been used to compare the values of 
(AA/A) 2 at the two frequencies. Booker [1958] lias 
pointed out that, in the Fraunhofer region (i.e., in 
the far-zone region, where (AA/A) 2 is independent of 
the distance between the observer and the diffracting 
screen), (AA/A) 2 should vary as the square of the radio- 
frequency, while in the Fresnel region (i.e., the near- 
zone region) the variation should be as the fourth 
power of the frequency. Figure 1 1 seems to indicate 
that, for values of (AA/A) 2 less than about 0.2 at 53 
Mc s, we are on the border between these two regions. 
For large values (not shown in the fig.) (AxA/A) 2 ap- 
proaches the Rayleigh limit (0.274) independent of 
frequency [Lawrence and Jespersen, 1959]. This is 
expected, since the Rayleigh limit represents the 
case where all of the signal (undiffracted component) 
has been removed. 

As has been indicated, the noise level may be esti- 
mated from the power spectra. This estimated Level 
has been used to compute the contribution of noise to 
t he value of (AA/A) 2 as a fund ion of (AA/A) 2 at t he two 
frequencies. The observed values of (AA/A) 2 were 
then corrected for noise effects and these are the 
values plotted in figure 1 1 . 

A rough indication of the diurnal variation, mixed 
with an unknown amount of seasonal variation, of 
amplitude scintillations comes directly from the sub- 
jective depth indices. The lower curve in figure 7 
was drawn by smoothing the hourly average values 
of depth index, and converting the result to (AA/A) 2 
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Figure 10. The Rice distribution function plotted in terms of 
power. 
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Figure 11. Amplitude scintillations at the two radiofrequencies. 

The expected wavelength dependence for Fraunhofer (X 2 ) and Fresnel (X 4 ) 
diffraction is shown hy the solid lines. 

by means of the table (see appendix 1). These 
curves are certainly too low because of a limitation 
at the high end of the scale of depth index. Index 5 
has been used to describe records in which the scin- 
tillations were "so deep that the lobes disappeared". 
Since many different degrees of scintillation activity 
are included in this category, the effects of the most 
violent scintillations are not properly represented 
in the low^er curves of figure 7. In order to assess 
the importance of this saturation effect, a portion of 
the data has been analyzed by the more elaborate 
procedure of evaluating P n (see following paper). 
The results, when transformed to equivalent val- 
ues of (AA/A) 2 and then averaged by months, 
are shown by the upper curves in figure 7. These 
curves give a more accurate estimate of the true 
variation of (AA/A) 2 . The data used in this result 
were a composite of two groups centered around the 
elevation angles 15° and 35°. 

Figure 12 shows the elevation-angle dependence of 
(AA/A) 2 as computed from P n . Only data from 
August 1958 to February 1959 are involved here. 
For comparison, a curve proportional to sec 3 x, the 
cube of the secant of the zenith angle, is included on 
the same figure. The observed variation does not 
agree very well with the sec 3 % curve which is pre- 
dicted by a theory based upon plane-earth geometry 
and spherical irregularities [Booker, 1958]. How- 
ever, most observers agree that the irregularities are 
alined along the earth's magnetic field and are 
there 4 fore not spherical. Also at the elevations con- 
sidered, the curvature of the earth is not negligible. 
Perhaps the most important reason for the dis- 
crepancy is that (AA/A) 2 approaches some constant 
value independent of frequency, and hence, for 
large scintillations, any theoretical elevation-angle 
dependence should be modified to include this effect. 

Using the noise-to-signal ratio, k, which indicates 
the fraction of incident energy that suffered diffrac- 
tion in the ionosphere, we can make a rough compari- 
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ELEVATION ANGLE (ABOVE NORTHEAST HORIZON) , DEGREES 

Figure 12. The variation, with source elevation above the 
northeastern horizon, of the amplitude index for each radio- 
frequency. 

The broken lines show a theoretical elevation angle dependence proportional 
to the cube of the secant of the zenith angle. 

son of our results with those of other observers. 
Interpolating between our two frequencies, we find 
that, for a given elevation angle, the value of k ob- 
served at Boulder is 2 or 3 times greater than the 
value observed at 80 Mc/s at Cambridge, England 
[He wish, 1952]. Since measurements at College, 
Alaska [Little et aL, 1958] have shown little or no 
dependence upon zenith angle, direct comparison 
with the Boulder results is less meaningful. The 
value of k at College was consistently 2 or 3 times 
greater than that at Boulder, even near the horizon, 
despite the fact that the College measurements were 
made at 223 Mc/s. (The Boulder and College ob- 
servations were made when the relative sunspot 
number was about 190; for the Cambridge observa- 
tions it was 85.) 

Angular scintillations. When normalized to ac- 
count for the spacing of the interferometer lobes, the 
phase plots can be used to calculate power spectra 
of angular scintillations over the frequency range 
0.005 to 0.17 c/s. Figure 13 shows examples of such 
spectra. The ordinate value is the common loga- 
rithm of power density, or variance per unit fre- 
quency interval, in units of tenths of a minute of 
arc, squared. Here, as with the corresponding 
spectra of amplitude scintillations, the reader is re- 
minded that spectra can be obtained only for moder- 
ately disturbed records. Records containing intense 
scintillation activity cannot be used to derive spectra. 

He wish [1951] has shown that for meter wave- 
lengths, the scale of the ionospheric irregularities is 
the same as the scale of the phase irregularities on the 
ground whenever the rms phase variations are less 
than 1 radian. That this is generally true for our 
data is shown by the fact that the shape of the power 
spectra of the angular scintillations are similar for 
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FREQUENCY, c/s 

Figure 13. The highest and I owes/ observed angular power 
Spectra for each radiof requeue if. 

The lower level at each frequency indicates the background "noise" level. W p 
is the variance per unit frequency Interval expressed in tenths of a minute of arc, 

squared. 2 is indicated, in units of minutes of arc squared. 



records taken simultaneously at the two radiofre- 
quencies. However, as in the case of amplitude 
scintillations, the high-frequency end of the power 
spectra is at times masked by the presence of noise 
which is folded hack into the spectrum. Figure L3 
shows the highest and lowest observed levels of the 
power spectra for angular scintillations at the two 
frequencies. As for amplitude scintillations, the 
lower two curves are indicative of the noise level. 
The corresponding values of the variance of the 
angular scintillations, # J , are given for each of the 
spectra. 

The shape of the spectrum varies sufficiently from 
day to day to obscure any variation in spectral shape 
which may be related to the elevation angle of the 
source. Inspection of simultaneous power spectra of 
angular scintillations at 53 and 108 Mc/s indicates 
that the details of the spectra do not coincide. The 
general shape may vary from hour to hour, usually 
in the same way on both radiofrequencies. 

Inspection of the phase plots of records taken 
simultaneously at the two radiofrequencies indicates 
(except for a special case to be described) that there 
is no detailed correlation. Here again, where dif- 
fraction is involved, no correlation is expected. 

Forty-one pairs of simultaneous phase plots (repre- 
senting moderate scintillation activity) were used to 
compare 2 , the mean square angular scintillation at 
the two frequencies. Figure 14 displays this fre- 
quency dependence of angular scintillations. 




Figure 14. Mean square angular scintillation at the two 
radiofrequencies. 
The solid line is the theoretical wavelength dependence. 
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ELEVATION ANGLE (ABOVE NORTHEAST HORIZON) , DEGREES 

Figure 15. The variation, with source elevation above the 
northeastern horizon, of the rms angular scintillation for each 
radiofrequency . 

The broken ijn rs show ,i theoretical elevation-angle dependence proportional 
to the square root of the secant of the zenith angle. 



Theory [Booker, L958] indicates that (here should 
be a fourth-power wavelength dependence for d 2 , and 
our data (its this quite well, Here there is no tend- 
ency to approach a limiting value which is independ- 
ent of frequency, since 6 2 depends upon the number 
of ionospheric irregularities which the radio wave 
passes through, each irregularity introducing a cer- 
tain phase change dependent upon the frequency of 
the incident wave. 

Figure 15 shows the elevation-angle dependence of 
6 2 at the two frequencies and figure 16 gives the 
seasonal (and/or diurnal) dependence of d 2 over a 
limited portion of the year. Both of these curves 
are based upon a (l (see following paper by the same 
authors) which was then transformed to the appro- 
priate value of 6 2 Simple theory indicates that d 2 
should vary as sec 1/2 x and this seems to be true in 
figure 15. 

Relationships between amplitude and angular scintil- 
lations. About two hundred scatter diagrams of the 
instantaneous power, P, versus the instantaneous 
phase, a, have been examined. In no instance is 
there an indication of a detailed correlation between 
these two quantities. This is the result to be ex- 
pected if a diffraction process is involved. However, 
a definite correlation does appear for the magnitude 
of power and phase variations averaged over 10 min 
intervals. 
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Figure 16. The seasonal and for diurnal variation of rms 
angular scintillation. 



Theory [Booker, 1958] indicates that the statistical 
relation between amplitude and angular scintillations 
depends upon the observer's distance from the iono- 
spheric disturbing region. In the Fraunhofer region 



whereas in the Fresnel region 



V A ) 8(xZ) ! 
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(1) 



(2) 



L is the ionospheric irregularity size and z is the 
observer's distance from the disturbing region . These 
results are valid only for rms phase scintillations less 
than 1 radian. From (1) it is evident that the rela- 
tion between (AA/A) 2 and 6 2 is linear, with slope pro- 
portional to (L/\) 2 in the Fraunhofer region. In the 
Fresnel region the relation depends upon z, the ob- 
server's distance from the disturbing region, and, 
furthermore, a plot of (AA/A) 2 versus d 2 should coin- 
cide for two different radiofrequencies. 

Figure 17 shows the relation which was obtained 
by fitting a line through the observed scatter of 
points. Because (AA/A) 2 increases with zenith 
angle, the horizontal scale corresponds to a non- 
linear scale of increasing z. The curves diverge, 
indicating a dependence on X, but they are also 
nonlinear, indicating a dependence on z. We must 
conclude that our observations were not completely 
within either the Fraunhofer or the Fresnel regions, 
but in the transition zone. 
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Figure 17. The observed relationship between mean square 
angular scintillation and amplitude index. 



Examination of the spectra of the power and phase 
scintillations from a given record indicates that there 
is no detailed correlation between the two, but that 
they have the same general shape, viz, steep at the 
low-frequency end, and leveling off toward higher 
frequencies. This is in agreement with other 
observations that the rate of fading of the amplitude 
and phase are similar for a given record [Ryle and 
Hewish, 1950]. 

Slow, irregular angular fluctuations. In addition 
to the angular scintillations described above, slow 
irregular phase variations having autocorrelation 
times of some minutes have been observed at both 
frequencies. Figure 18 shows two typical examples 
of these irregularities as observed simultaneously at 
the two radiofrequencies. If one were to remove 
the slow variation from them they would look like 
ordinary scintillation records in which no irregu- 
larities are present. Figure 19 shows one of the slow 
irregular structures measured for a period of 1 hr. 
This figure is a composite of consecutive records 
of the type shown in figure 18. It should be noted 
that the steady angular variation of the radio star, 
due to its diurnal motion, has been removed, as it 
is from all phase plots. Since no attempt has been 
made to determine the absolute phase of the inter- 
ferometer system, the zero levels are arbitrary on 
these diagrams. 

The slow angular variations show high correlation 
between records taken simultaneously at the two 
frequencies. Forty-seven pairs of records, each ten 
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Figure 18. Typical examples of slow irregular angular varia- 
tion observed simultaneously at the two radiofrequencies. 
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Figure 19. A sixty minute record exhibiting slow irregular 
angular variation. 
108 Mc/s; 1100 a.m., October 17, 1958. 



minutes long, showing slow irregular variations at 
the two frequencies, have been examined to deter- 
mine the average total angular excursion for each 
frequency. The ratio of the average excursion at 
53 Mc/s to that at 108 Mc/s ranges from 2:1 to 6:1 
with the average 4.1:1. A four-to-one ratio would 
be expected if ionospheric refraction were responsible 
for the angular excursions. The close agreement 
between the observed ratio and the theoretical 
ratio expected for refraction, together with the high 
cross-correlation, indicates thai these slow angular 
variations are due to ionospheric refraction. Tests 
with a target transmitter were made to eliminate 
the possibility that equipment malfunction might be 
responsible for these slow angular variations. The 
negative results of these tests, together with the 
frequency dependence of the variations, indicate 
that the effects are of ionospheric origin. 

Examination of the amplitude variations of pairs 
of records for which there were slow angular scintilla- 
tions, indicates that there is no correlation between 
the amplitude fluctuations on the two frequencies, 
nor is there any long-period variation in amplitude. 
Jn fact, the amplitude scintillations seem to be un- 
affected by the presence of the slow angular 
scintillations. 

This lack of correlation suggests that the iono- 
spheric irregularities responsible for the refraction 
are acting as prisms, or as weak lenses whose focal 
lengths are long compared with the observer's 
distance from them. For the case of a prism we can 
use the observed total angular excursion to estimate 
the amount by which the electron-content of the 
ionosphere deviates from its mean value. Smith 
[1952] has shown thai a prism will introduce an 
excursion in apparent direction given by 

e 2 dn 
Z 2irmP dx 

Here dn/dx is the gradient, transverse to the line of 
sight, in electron content per unit column along the 
line of sight; /is the radiofrequency ; e is the charge 
of an electron ; and m is the mass of an electron. 

A typical angular excursion at 53 Mc/s is about 
0.5 (leg. If we assume that such a gradient is 
maintained over a portion of the irregularity equal 
to, say, 10 km, then from the above formula, the 
total electron-content variation is a few percent of 
the mean value. 

Figure 20 shows the seasonal and or diurnal varia- 
tion in the frequency of occurrence of these irregu- 
larities. The ordinate value is the percent of obser- 
vations which exhibit slow irregular variations, 
relative to the total number of cases available for 
each month. The shape of the curve suggests that 
the variation is more diurnal than seasonal, since 
the peak is centered about two seasons (fall and 
winter). Furthermore, in the months from May 
to December the data were taken during daylight 
hours. The shape of the curve seems to indicate 
that the irregularities build up rapidly as the iono- 
sphere becomes illuminated, and then gradually 
disappear. 
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Figure 20. The seasonal and/or diurnal variation in the 
frequency of occurrence of slow irregular angular varia- 
tions. 



The total angular excursion is found to decrease 
as the source rises, while the probability of observing 
these slow irregularities seems to increase. Since 
ordinary amplitude and phase scintillations decrease 
as the source rises, the slow irregularities are most 
likely to be seen on records for which the ordinary 
scintillations are small, (AA/A) 2 less than 0.06. 
However, they have occurred on several occasions 
when (AA/A) 2 was as large as 0.20. 

It is reasonable to assume that the elevation- 
angle dependence in the probability of occurrence 
of the slow angular variations results from the fact 
that at low elevation angles amplitude scintillations 
are most severe and are likely to limit the ability of 
the " zero-crossing" method to recover the slow 
phase variations from the data. 

The variation in the total angular excursion with 
elevation angle probably arises from the changing 
path length through the ionosphere, although there 
is some possibility that the structure of the irregu- 
larities^jairies rapidly with latitude. No attempt 
has yet been made to check these ideas quantita- 
tively. 

At the low elevation angles of 10 to 20 deg, the 
distance which the line of sight moves through the 
ionosphere in a period of 20 min (typical of the length 
of time during which any one irregularity is observed) 
is of the order of 200 km. This distance may be 
assumed to be indicative of the horizontal size of the 
ionospheric irregularities in the appropriate direc- 
tion if ionospheric drifts are neglected. 

Examination of the distributions of ionospheric 
parameters scaled from the Ellsworth ionograms 
reveals no significant difference between times of 
slow angular scintillations and other times. A 
possible exception is a marginally significant increase 
mfoF2. 



Vitkevich and Kokurin [1957 and 1958], using a radio- 
astronomy technique, observed large-scale angular 
variations with no corresponding amplitude effects. 
Other authors have apparently observed a similar 
phenomenon [Bramley and Ross, 1951; Whitehead, 
1956; Jones et al., 1957]. The large-scale travel- 
ing disturbances which Munro [1950] has studied are 
possibly related to the disturbances observed here. 
Satellite observations made at the Boulder Labora- 
tories of the National Bureau of Standards [Little 
and Lawrence, 1960] have also disclosed the presence 
of these large-scale, ionospheric irregularities. 

6. Conclusions 

After removal of the diurnal (and/or seasonal) 
trends from the data, amplitude scintillations still 
correlate with spread F. The variation of correla- 
tion with elevation angle of the source is small, 
indicating that spread F is in patches whose sizes 
are at least 400 km. No significant correlation is 
found between amplitude scintillations and sporadic 
E. This result is particularly interesting since the 
experiment was designed especially to measure and 
detect possible E'-region effects on scintillations. 

The distribution of flux received from a scintil- 
lating radio source closely matches the theoretical 
distribution to be expected from the superposition of 
several randomly phased, diffracted waves upon a 
constant, undiffracted wave. The zenith-angle de- 
pendence of these amplitude scintillations does not 
agree with a theory based upon isotropic ionospheric 
inhomogeneities. At Boulder, the fraction of the 
radio-star energy diffracted by the ionosphere is 
about twice as great as at Cambridge, England, but 
less than half as great as at College, Alaska. How- 
ever, the Boulder and College observations were 
made when the relative sunspot number was about 
190; for the Cambridge observations it was 85. 

The root-mean-square value of angular scintil- 
lations is found to be proportional to the square of the 
wavelength, in accord with a theory of diffraction by 
ionosph eric irregularities. 

Comparison of amplitude scintillations with angu- 
lar scintillations indicates that, for elevation angles 
of 15° to 50°, the region of the ionosphere responsible 
for scintillations lies near the border between the 
Fresnel and Fraunhofer diffraction regions for both 
frequencies. 

Slow irregular angular variations having autocor- 
relation times of 5 to 20 min are commonly observed 
at both frequencies. These variations occur chiefly 
by day and are attributed to large (of the order of 
200 km) lens-like ionospheric irregularities whose 
electron content per unit vertical column varies by a 
few percent from the mean value. The resulting 
angular disrla cement of a radio source is frequently 
0.5 deg at £3 Mc/s, and varies as the square of the 
wavelength. 

7. Appendix 1 

Amplitude scintillation index. A visual estimate of 
the depth of scintillations can be made from a pen- 
recorder trace of the interferometer output. Figure 
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21 shows L2-min sections of two typical scintillation 
records. The modulated sine wave is formed as the 
interferometer lobes sweep past the discrete source 
(Cygnus-A) at the rate of 10/min. The amplitude 
of the sinusoidal trace is proportional to the flux 
received from the discrete source; the fluctuations in 
the envelope are ionospheric scintillations. The 
depth of these amplitude scintillations was classified 
by eye into 5 categories as follows : 



Depth 
index 



Appearance of record 



No scintillations. 
Shallow scintillations. 
Medium scintillations. 
Deep scintillations. 
Scintillations so deep that 
the lobes disappear. 



For example, the records in figures 21a and 21b have 
been assigned depth indices 3 and 5 respectively. 



The estimated depth index for each 5-min period was 
inserted upon the appropriate IBM card. 

Many of the scintillation records were supple- 
mented by simultaneous digital recordings on 
punched paper tape, A random sample of 108 Mc/s 
digital records was selected for each depth index and 
subjected to detailed analysis on an electronic com- 
puter. This analysis yielded both the power index 



AP = 

P = 

and the amplitude index 



P-P 



(^)W)* 



where P is the power received from the discrete 
source, A = ^P is the amplitude of this radio signal, 
and a bar denotes an average over the record, 
results are summarized in the following table: 



The 







AP/P 






(AA/A)* 




Depth 














index 
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nun 
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Figure 21. Twelve-minute samples of typical 108 Mc/s phase-sweeping records. 

The modulated sine wave is formed as the interferometer lobes sweep past the source. The amplitude of the record is proportional to the flux received from the- 
source; the fluctuations arc ionospheric scintillations. 
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8. Appendix 2 

Correlation of amplitude indices with ionograms. 
The IBM cards containing the simultaneous iono- 
spheric (Ellsworth) and scintillation data are referred 
to as "statistical index cards". Tests for correlation 
between any two parameters, x and y, on the cards 
may he performed easily on IBM statistical machines. 
The statistical-sorter output consists of a bivariate 
distribution table, each variable, x and y, being 
allowed a scale of one through nine, zero being re- 
served for cases where there are no data. These 
ten-by-ten distribution arrays are then punched onto 
IBM cards for processing through a computer pro- 
gram which determines the correlation coefficient 
between the two variables. The normalized corre- 
lation coefficient is defined in the usual manner as 



r= 



xy—xy 



(i) 



where the variance <t 2 x =x 2 — x 2 . For independent 
data which are normally distributed, the standard 
deviation of the correlation coefficient is 



1-r 2 



(2) 



where N is the number of samples. If the data are 
not independent, but possess a significant serial 
correlation over an interval n samples in length, then 
iVin (2) should be replaced by N/n. A visual inspec- 
tion of the data indicates n to be approximately 6 
for data taken at 5 min intervals, but nearly 1 for 
data taken at half-hour intervals. For the correla- 
tion coefficients reported here, the standard devia- 
tions have been calculated using these values of n. 

The distribution of the data is approximately 
normal for most of the ionospheric parameters but 
for scintillations it is strongly skewed towards the 
high depth indices. Figure 22 gives the distribution 
of scintillation index for each radiofrequency, both 
before and after revision of the data. The methods 
of, and reasons for, data revision will be discussed 
in the following paragraphs. In order to assess the 
error introduced by using (2) with skewed data, a 
more sophisticated analysis [Cramer, 1954] was 
undertaken. This calculation, involving higher order 
moments of the actual distribution, was made on a 
typical correlation coefficient taken from our data. 
The difference in o> amounted to less than 0.001; 
hence (2) has been used in all subsequent calculations 
of a r . 

In cases where \r\ >0.5, or where iVis small, there 
is some advantage to use of the Fisher Z r function 
[Walker and Lev, 1953] and its corresponding stand- 
ard deviation a z . Z r is defined as 



z =l ln Q^)' 



(3) 
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Figure 22. Frequency distribution of observe^ scintillation 
index. 



and cr, is obtained from the formula 



v W-3 



(4) 



In the present case where correlation is small and 
N is large, Z r and a Zr are practically equal to r and 
o> respectively, and they can be used interchangeably. 

The geometry of the experiment indicates that 
any correlation between ionospheric and scintillation 
parameters should vary with local sidereal time 
(LST) as the ray path from Cygnus-A moves through 
the different effective ionospheric levels over Ells- 
worth. Also, scintillation activity is found to 
decrease with increasing elevation angle (increasing 
LST) of the source. Consequently, the statistical 
index cards were divided into groups according to 
the elevation angle of Cygnus-A. Figure 3 gives 
these groupings as well as a complete distribution 
table of the cards as a function of the elevation angle 
of Cygnus-A. Figure 3 also gives the relationship 
between the elevation angle of Cygnus-A, local side- 
real time (LST), and effective ionospheric heights. 
The average elevation angle of each major (10°) 
group has been computed and is used for plotting 
the correlation coefficients in section 5. 

Early correlation studies, reported previously 
[Lawrence and Jespersen, 1959], consisted primarily 
of bivariate correlations as a function of LST, and 
used only six months' data. With the cessation of 
recording, the thirteen months (February 1958 
through February 1959) of available data were 
placed on statistical index cards. These data were 
then tested for correlations, but with few consistent 
results. By contrast, the previously reported 
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results [Lawrence and Jespersen, 1959] of correlation 

of a portion of the data had appeared internally 
consistent. Prominent diurnal and/or seasonal 
effects in the full year's data are presumably 
responsible for this difference. 

Division of the data into two classes, sunlit and 
dark ionosphere, proved helpful and indicated that 
complete removal of diurnal and seasonal trends 
would be worthwhile. To determine the trends, the 
hourly average value of each variable w T as plotted, 
except that in the case of spread F and sporadic E 
the frequency of occurrence was used, as shown in 
figure 6. The hourly average values of scintillation 
depth index, converted to (AA/A) 2 , appear in figure 7. 
These graphs are referred to as diurnal curves, 
although each contains an unknown amount of 
seasonal variation. It is impossible to distinguish 
with certainty between diurnal and seasonal effects, 
since only one radio source was used. 

The IBM ()50 computer was used to subtract the 
appropriate hourly average 4 from each numerical 
value on the statistical index cards, and to normalize 
the remainder to a scale of one through nine. The 
resulting cards are referred to as "revised statistical 
index cards", and have been used in all the correla- 
tion results reported here. Figure 22 shows the 
effect of this revision upon the distribution of the 
scintillation indices. This redistribution cannot, of 
course, remove the effect of skewness in the original 
depth indices, as this was a result of limited equip- 
ment sensitivity. 

The data in final form consist of several thousand 
revised statistical index cards. The extreme excur- 
sions of the inospheric and scintillation indices from 
their average value for any particular local time, 
now appear as extreme values on a scale of one 
through aine. Most of the data are clustered about 
the median of the scale and so they contribute little 
to the correlation coefficient. Thus the revised 
values are in a form which may be tested for detailed 
correlation between individual events, regardless of 
trends in the averages of the original variables. 



This experiment could not have been performed 
without the cooperation of several people. Among 
them are E. R. Schiffmacher and H. A. Erickson 
who were responsible for the maintenance and 
operation of the scintillation receivers, It. F. Carle 
who established and operated the ionosonde at 
Ellsworth, and Mrs. C. I. Enfield who assisted w T ith 
the analysis of t he data. This project was supported 
by the Ballistic Missiles Division of the United 
States Air Force under D. O. AF 04(647)-134. 
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